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ABSTRACT 

The physical mechanisms that set the initial rotation rates in massive stars are a crucial unknown in 
current star formation theory. Observations of young, massive stars provide evidence that they form 
in a similar fashion to their low-mass counterparts. The magnetic coupling between a star and its 
accretion disk may be sufficient to spin down low-mass pre-main sequence (PMS) stars to well below 
breakup at the end stage of their formation when the accretion rate is low. However, we show that 
these magnetic torques are insufficient to spin down massive PMS stars due to their short formation 
times and high accretion rates. We develop a model for the angular momentum evolution of stars over 
a wide range in mass, considering both magnetic and gravitational torques. We find that magnetic 
torques are unable to spin down either low or high mass stars during the main accretion phase, and 
that massive stars cannot be spun down significantly by magnetic torques during the end stage of 
their formation either. Spin-down occurs only if massive stars' disk lifetimes are substantially longer 
or their magnetic fields are much stronger than current observations suggest. 

Subject headings: stars: formation - stars: magnetic field - stars: massive - stars: protostars - stars: 
rotation 



1. INTRODUCTION 

While there has been significant theoretical attention 
to understanding the initial rotation rates of Sun-like 
stars, far less work has been done on more massive 
stars. Since the stellar evolutionary path depends on 
the rate of mass loss and i nternal mixing, both of which 
are enhanced by rotatio n (|Biorkman fe CassinellilH993t 
iMaeder fc Mevnetll2010l) . our inability to predict initial 
rotation rates is a limiting factor in stellar evolution the- 
ory. Observations of young, massive stars provide evi- 
dence that they form in a similar fashion to their low- 
mass counterpa rts: via gravitational collapse of a molec- 
ular cloud co re (iMcKee fc Tanll200l iZapata et aI1l2008t 
iDavies et all 1201 ID . These cloud cores are slowly rotat- 
ing but have very large radii, and thus have high ini- 
tial angular momenta. This has led to the "angular 
momentum problem" in which the initial angular mo- 
mentum of a cloud core is at least thre e orders of mag- 
nitude greater than the resulting star (jGoodman et al.1 
Il993t lBodenheimerlll995t lLarsonll20l6T) and must be re- 
distributed or removed during collapse. 

Massive stars f orm in magneti zed high-density turbu- 
lent gas clumps (jCrutcherl 119991 ) that are characterized 
by short core collapse times and h igh time-averaged ac- 
cretion rates (McKee & Tan 2003). Due to the high an- 
gular momentum content of the diffuse gas, material is 
unable to be directly deposited on to the central object 
and is instead circu larized at a distance fa r from the star, 
resulting in a disk (jKrumholz et all [20071 l2009h . Obser- 
vations, although rare, confirm that dis ks form around 
massi ve pro tostars durin g clou d collapse (jCesaroni et al.l 
120061 120071: iChini et al.l l2011h and the accretion onto 
these disks is regulated at le ast in part by the magnetic 
field (|Vlemmings et ah! 1201 Of ). Furthermore, these disks 
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might evolve like those located around young, low-mass 
stars (|Chini et al.1 120061 ). The disk transfers mass and 
angular momentum to the central protostar, which acts 
to spin it up. This transfer of angular momentum, along 
with contraction of the protostar towards the main se- 
quence, suggests that young stars should be rotating at or 
near their break up speed, the rotational speed at which 
t he centripetal for ce at the equator balances gravity. 

iLin et aD ()2011l ) found that gravitational torques pro- 
hibit a star from rotating above ~ 50% of its break 
up speed during formation. However, the observed pro- 
jected rotation rates of young low mass and some mas- 
sive stars suggest that they rotate at a much lower frac- 
tion. Observations of low-mass PMS stars suggest that 
their rotation periods span a factor of ~ 30 and ap- 
proximately half are slow rotators, rotating at about 
10% of their break up speed (jHartmann fc Staufferlll989l : 
iHerbst et al.ll2007l) . The observed rotational velocities of 
massive stars suggest that they are spinning significantly 
faster than their low-mass counterparts. iWolff et al.l 
(|2006f ) studied a sample of young massive stars (M* > 
25 Mq) and found that their median rotation rate was 
20% of their break up speed. iHuang et al.l (|2010D ob- 
served the projected rotational velocity distribution of 
220 young B stars and found that approximately 53.3% 
are rapid rotators, rotating with a velocity that is at least 
40% of their break up speed. How these initial rotation 
rates are achieved and their dependence on stellar mass 
is still an unanswered question. 

The physical mechanism responsible for causing young 
low-mass stars to be slow rotators has received consid- 
erable attention over the last three decades. One pop- 
ular theory is that during the T Tauri phase (experi- 
enced by PMS stars with masses less than ~ 3 M Q ), 
when the accretion rat e is low, M a < 10 -7 M Q yr _1 
(jHartmann et all I2006D . the magnetic connection be- 
tween the star and its accretion disk can transport sub- 
stantial angular momentum away from the star, result- 
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ing in spin rates well below bre a k up in agreement 
with observations (|Koenigll 119911 : lArmitage fc Clarke] 
ll996j). The fact that T Tauri stars have strong mag- 
neti c fields, typically, several hundred G to several 
kG (j Johns-Krulll 120071 ) . long contraction timescales af- 
ter their main assembly, and long acc retion disk life - 
times su pport this spin down s cenario (iBouvier l 12001 . 
However, IMatt fe Pudritd (|2005t) and iMatt et al.l (|2010f ) 
found that when the stellar magnetic field lines open due 
to the differential twisting between the star and disk the 
resulting rotation rates, while still below break up, are 
higher than those of the slowest rotators. 

Magnetic fields have been detected in a small sample 
of young and evolved OB stars. These fields arc between 
a few hundred G to several kG and typically have a 
bipolar topology (iDonati et al.ll200l iWade et all 120061: 
Hubrig et all 120081 : iGrunhut et al.1 120091 : iMartins et al.1 



20101) . The origin of these fields is poorly understood, 



since the envelopes of such stars are radiative rather 
than convective, ex cluding the possibility of a Solar- 
type dynamo effect (|Mosa 12001). The favored hypoth- 
esis for the presence of magnetic fields in massive stars 
is that they are fossil fields that wer e either accumulate d 
or generated during star formation (| Warder et al.ll20lTI ). 
lAlecian et al.l (|2008l ) discovered two very young B stars 
with strong surface magnetic fields. They found that the 
younger of the two is a rapid rotator and situated in the 
first half of the PMS phase, whereas the older star, which 
might already be on the main sequence, is a slow rotator 
most likely spun down via magnetic torques. 

This implies that massive stars likely have strong mag- 
netic fields present during their formation and that these 
fields, due to coupling with the accretion disk, may be 
able to remove a substantial amount of angular momen- 
tum from the star, producing spin rates on the zero- 
age main sequence (ZAMS) well below break up in a 
similar fashion to their low-mass counterparts. How- 
ever, massive stars reach the ZAMS very quickly since 
they have short thermal equilibrium timescales. They 
also have higher accretion rates during their formation 
and their magnetic fields are weaker relative to their 
stellar binding energy as compared to low mass stars. 
They likely have shorter disk lifetimes than contract- 
ing low mass stars, since their disks are likely to be 
quickly photo-disinteg rated due to their high luminosities 
(jCesaroni et al.l 120071 ) . All of these factors make mag- 
netic spin-down more difficult. In this paper we explore 
whether the initial spins of massive stars are regulated 
by the interaction of their accretion disk with the stellar 
magnetic field. To study this issue we model the angular 
momentum evolution for both low-mass and massive pro- 
tostars by considering both magnetic and gravitational 
torques. We apply the star-disk interaction model devel- 
oped bv IMatt fc Pudritd (|2005f ) (hereafter MP05) where 
the stellar magnetic field is connected to a finite region of 
the accretion disk, and the twisting of the magnetic field 
lines due to the differential rotation between the star and 
disk leads to a spin-down torque on the star. 

This paper is organized as follows. In the following sec- 
tion (§2), we give a brief introduction to how the presence 
of surface magnetic fields during the protostcllar phase 
can extract angular momentum from the star. We de- 
scribe our stellar angular momentum evolution model, 
which include a prescription for protostellar evolution 



and the star-disk interaction, in §3. We state our results 
in §4. Lastly, we discuss our results in §5. 

2. MAGNETIC TORQUES: THEORY & BACKGROUND 

Protostars embedded in circumstcllar disks accrete ma- 
terial from an angular momentum-rich mass reservoir. If 
the disk is Keplerian the specific angular momentum con- 
tent of the circulating material, j = y/GM+r, increases 
outward and the angular velocity increases inwards. The 
presence of a stellar magnetic field is able to disrupt the 
disk outside the stellar radius and channel the disk mate- 
rial along field lines. Spin-down torques will be conveyed 
to the star due to the differential twisting of the field 
lines threading the accretion disk at radii where the disk 
rotates at a lower rate than the star. In this section we 
give simple scaling arguments to demonstrate how spin 
evolution varies with stellar mass, before proceeding to 
a more detailed numerical model in §3. The derivation 
that follows is an oversimplification and ensures maxi- 
mum spin down via magnetic braking. We include this 
section for the reader who is unfamiliar with the litera- 
ture. 

The radial extent of the accretion disk can be altered 
if the protostar has a magnetic field. The magnetic field 
is able to truncate the disk at the Alfven radius (denoted 
Ra) where the magnetic pressure, B 2 /8tt, balances the 
ram pressure, pv 2 , of the infalling material. Assuming 
the stellar magnetic field is dipolar and the magnetic 
field axis is aligned with the rotation axis of the star, 
the z component of the field in the equatorial plane at a 
distance r from the star is given by 



B z — 



-3 



(1) 



where B+ is the magnetic field strength at the stellar 
surface. The location at which the magnetic pressure 
is able to truncate the disk, assuming spherical free-fall 
accretion, is 



10- 7 M Q yr" 1 
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where M a is the accretion rate. In the case of disk accre- 
tion, the truncation radius is in general smaller than the 
value given in equation (2) by a factor of order unity. For 
simplicity and for the purpose of this section we neglect 
this factor in the following discussion. 

If the stellar magnetic field lines are connected to the 
disk the differential rotation between the two will cause 
the field lines to twist in the azimuthal direction induc- 
ing torques on the star. The disk co-rotates with the star 

at the location R co = (GM*) 1 / 3 S\ 2 ^ 3 where is the 
angular velocity of the star. The stellar field lines that 
connect to the disk outside the R co spin up the disk and 
spin down the star. If the field lines connect to a signif- 
icant portion of the disk outside of R co the star can be 
spun down to a velocity well below its break up speed. 

The stellar magnetic field lines threading an annulus 
of the accretion disk with width dr will exert a torque: 



B*B*r*dr 



(3) 
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where B^ is the azimuthal component of the field gener- 
ated by the twisting of the field lines relative to the star 
and is given by 



B 



n(r) - n+ 

fi(r) 



(4) 



where Q is the angular velocity of the Kcplerian accretion 
disk. Integrating equation (3) from Ra to infinity the 
total torque on the star due to the stellar magnetic field 
lines connected to the disk is 



{ * K * ( P-3 9R -3/2 R -3/2 

— — yu A - ZH co H A 



(5) 



The accretion of disk material at Ra adds angular mo- 
mentum to the star at a rate 



= M a ^GM*R A . 



(6) 



Notice that equation (5) contains both spin-up and spin- 
down torques acting on the star due to field lines con- 
nected to the disk within and outside of i? co , respectively. 
In order for the net magnetic torque to transport angular 
momentum away from the star (i.e., r m < 0) Ra must 
be greater than 



Ra, 



0.63i? c 
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In a system where the stellar parameters (M*, i?*, 
B±, M a ) are relatively constant there exists an equi- 
librium state, called the "disk-locked" state jKoenigll 
flflMtlArmitage fc ClarkellQQfitlMatt fc Pudritzll2005h . in 
which the stellar spin rate will adjust to its equilibrium 
value (i.e., when r a + r m = 0). Setting r a = — T m the 
equilibrium spin rate, as a fraction of the break-up speed 

(^bu 
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Assuming that the moment of inertia of the star stays 
constant, the characteristic timescale to reach equilib- 
rium is: 



U, eq = k 2 M*Rl 



(9) 



where k is the dimensionless radius of gyration whose 
value depends on the stellar structure. Equation (7) only 
holds when Ra > R*, which is true if the star has a 
surface magnetic field strength above a minimum value: 



> 400 



M a 



1/2 

IQ- J Mg yi— 1 / 
1/4 / n x -5/4 
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Figure [T] shows the equilibrium spin rate as a frac- 
tion of the star's break up speed and the correspond- 
ing time scales required for a 1 M Q star and a 30 
M Q star to reach equilibrium starting from rotation at 
break up, both as a function of the accretion rate. We 
adopt surface mag netic field streng t hs of 2 kG simi- 
lar to observations (|Wade et al.ll2006l: Uohns-Krulill2007l: 
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Fig. 1. — The equilibrium spin rate of a star as a fraction of its 
break up spin rate (top panel), and the corresponding spin down 
time scale (bottom panel) for 1 Mq (black solid line) and 30 Mq 
(teal dashed line) stars to reach equilibrium. Both stars have a 
surface magnetic field strength of 2 kG with a dipolar topology 
and are initially rotating at break up. The horizontal line in the 
top panel shows where the equilibrium spin rate is equal to the 
break up rate. 



iGrunhut et alj|2~009[) and assume k = 0.27 for a radiative 
star (e.g., n=3 polytrope). We adopt radii of 3 R Q for the 
1 Mq star (the typical radius of a T Tauri star of this 
mass) and 7.76 R Q for the 30 M Q star (ZAMS value). 
We consider only accretion rates where the equilibrium 
spin rate is below the break up rate. As the accretion 
rate increases, the equilibrium spin rate approaches the 
break up rate and the equilibrium timescale quickly de- 
creases. We find that magnetic torques produce equilib- 
rium spin rates below break up only for accretion rates 
below M a < 5 x 10 -5 M Q yr _1 , regardless of the stellar 
mass. In this regard, low- and high-mass stars are sim- 
ilar. The typical mass accretion rates during the main 
accretion phase, where the majority of the stellar mass 
is accreted, for low- and high -mass star formation are 
5 x 10~ 6 vr" 1 jShul [l977h and 5 x 10~ 4 M Q yr" 1 
(|McKee fc Tanll2003D . respectively. For our adopted field 
strength, Ra for the 30 M Q star is within the stellar sur- 
face at this accretion rate. In contrast, the disk is trun- 
cated very close to the stellar surface for the 1 M Q star, 
leading to an equilibrium spin rate close to break up. We 
conclude that disk truncation does not occur for massive 
stars and is unimportant for low-mass stars during the 
main accretion phase. At the lower accretion rates that 
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are likely to occur after the main accretion phase ends, 
we find that low- and high-mass stars differ in that the 
latter have much longer equilibration timescales than the 
former due to their larger inertia. For example, the equi- 
libration timescale for the 30 M Q star for very low accre- 
tion rates is a significant fraction of i ts stellar lifetime, 
t ms = 5.9 Myr (|Parravano et al.ll2003| ). Furthermore, at 
high accretion rates this ti mescale is comparab le to the 
star's formation timescale ([McKee k Tadl2"o"03[) suggest- 
ing that massive stars are unable to reach spin equilib- 
rium. To further explore the consequences of this analy- 
sis we follow the angular momentum evolution of massive 
protostars to determine the physical conditions that are 
required to spin them down by magnetic torques. 

3. STELLAR ANGULAR MOMENTUM EVOLUTION 
MODEL 

The goal of this work is to determine if the initial rota- 
tion rates of massive stars can be regulated by magnetic 
torques due to the interaction of the stellar magnetic 
field and surrounding accretion disk during formation. 
To this end, we construct a simple model to track the 
mass, radius, and angular momentum content of accret- 
ing protostars subjected to gravitational and magnetic 
torques. We describe the elements of this model in the 
following subsections. 



torque on the star due to the coupling of the stellar mag- 
netic field with the surrounding accretion disk described 
in §3.3. 

3.2. Accretion History 

The accretion history of our protostars is divided into 
two distinct accretion phases. The first is the main accre- 
tion phase given by the turbulent core model from MT03, 
which describes an accelerating accretion rate, where the 
majority of the stellar mass is accreted. This model as- 
sumes that the star-forming core is marginally unstable, 
massive, and supported by turbulent motions. Next, we 
follow the disk clearing phase in which the accretion disk 
is no longer being fed by the core envelope. These accre- 
tion phases are described in §3.2.1 and §3.2.2. 

3.2.1. Primary Accretion Phase: Core Collapse 

We model the mass accretion using the two-component 
core model of MT03 which assumes the central region 
of a molecular cloud core is dominated by thermal mo- 
tions and the core envelope is d ominated by non-thermal 
motions (jMvers k FuHerlll992l ). This leads to a density 
distribution that is equivalent to the sum of a singular 
poly tropic sphere and a singular isothermal sphere: 



3.1. Protostellar Model 

We monitor the spin and angular momentum evolution 
by following the protostellar radius and internal structure 
evolution dur ing its formation wit h the use of the one- 
zone model of iMcKee k Tad (120031) (hereafter MT03) as 
updated by lOffner et al.l ( 2009(1 . By treating the proto- 
star as an accreting polytrope and requiring conservation 
of energy, the evolution of the protostellar radius is given 
by: 



2M a R„ 
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where M a is the accretion rate onto the protostar, /k 
is the fraction of kinetic energy of the infalling ma- 
terial that is radiated away, ftp is the ratio of radia- 
tion pressure to the total pressure, a g = 3/ (5 — n) is 
the coefficient describing the binding energy of a poly- 
trope, Li n t is the internal stellar luminosity, L\ is the 
rate of energy required to dissociate and ionize the in- 
falling material, and Lp> is the rat e at which energy i s 
supplied from burning deuterium (|Nakano et al.l I2000D . 
The model also includes a few discontinuous changes in 
polytropic index and radius to represent events such as 
the onset and cessation of core deuterium burning and 
the formation of a radiative core. We u s e the model 
parameters recommended by lOffner et al.l (|2009h which 
ar e based on the detailed stell ar evolution calculations 
by IHosokawa k Omukail (I2009D. We refer the reader to 
MT03 and Appendix B of lOffner etHI (pOOl for a de- 
tailed description of the model and protostellar evolu- 
tionary states. 

We treat the protostar as a solid body to follow its an- 
gular momentum content ( J* = 7*0 + ). We evolve the 
stellar angular momentum content by computing the net 
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where p s is the density at the surface of the core, i? cor c is 
the core radius, and cth is the thermal sound speed within 
the core and is assumed to be constant. We adopt the 
fiducial value of k p = 1.5 from MT03 in agreement with 
observations describing the turbulence-s upported density 
profile of massive star forming cores (|Caselli k Mversl 
[19951 Fvan~der Tak et al J 120001 iBeuTher et al.l^OOlt l 

The accretion rate onto the disk, which is supplied by 
the background core, is: 



<t>*M*,t 
t*.ft 
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(13) 
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where t^s = (37r/32Gp) is the free fall time evaluated 
at .Rcore, is the current stellar mass, M+.f is the final 
stellar mass, and 



3(2 -k p 
2(3- k 



(14) 



The dimensionless constants cj>, </>*.th, and </>*,nth are of 
order unity and depend on k p and the magnetic field 
strength. The efficiency factor, e coro , describes how much 
of the core mass will end up in the star rather than be- 
ing ejected by the protostellar outflow and we adopt the 
value of 0.5 from MT03, whic h is typical of both low- 
mass dMatzner k M cKcc 2 000l) a nd high-mass star for- 
mation ([Cunningham et al.ll201l[ ). The parameter Af tn 
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Fig. 2. — Accretion history of protostars with final masses of 
0.5 - 50 Mq, following equations (13) and (16), for our fiducial 
parameters given in table 1 in §4. 

describes the mass below which the thermal density dis- 
tribution dominates. For a core with surface density 
£ = M Jrif e~ 1 re /7ri?2 oro , M t h is defined as: 



1.23 x 10 



-3 ( T 1 
I 20 K) 



M, 



(15) 



where Sq = g cm -2 ). We further assume that the 
accretion rate onto the disk is the same as that onto the 
star and use this value for our protostellar accretion rate. 

3.2.2. Secondary Accretion Phase: Disk Clearing 

Late in the formation the core envelope will exhaust its 
reservoir of mass and no longer feed the accretion disk. 
We assume that we are left with a thin, Keplerian accre- 
tion disk that continues to transfer mass and angular mo- 
mentum to the central protostar. For simplicity and be- 
cause observations of disks located around massive stars 
are very limited, we assume that this results in a decreas- 
ing accretion rate as a fun ction of time which we model 
as a dec aying exponential (iCollier Cameron fe Campbell 
[1991 lYil fl99l [19951 iMatt et al.ll2010D ~ 



Mj 



(16) 



where Md is the remaining mass in the accretion disk 
(i.e., the total amount of mass that would accrete from 
t = — >• oo) and t a is the decay timescale. Since Md and 
tg, are highly unconstrained, we experiment with differ- 
ent values in §4. Figure [2] shows the accretion history, 
including both the core collapse and disk clearing accre- 
tion phases, for stars with final masses of 0.5-50 M©. 

3.3. Star-Disk Interaction Model 

In §2 we showed how the presence of a stellar mag- 
netic field can remove angular momentum from the star 
as it accretes matter from an accretion disk. This de- 
scription assumed that the stellar field lines were con- 
nected at all radii of the disk larger than Ra- However, 
the differential rotation between the star and disk will 



twist the connected field lines. This twisting will cause 
the magnetic field to undergo a rapid inflation leading 
to an opening of the field lines, effectively decreasing 
the size of the disk region that is connected to the stel- 
lar magnetic field (lLovelace et all 119951 : lUzdcnsk v et al.l 
I2002t IMatt fe Pudritzll2005l) . We now include this effect 
when calculating the net magnetic torque on the star 
with the use of the model developed by MP05, which 
is an extension to the dis k-locking model first developed 
bv lGhosh fe Lambl (119781) f or accreting neutron stars and 
extended bv lKoenigll (|199ll) to describe the star-disk cou- 
pling for magnetized T Tauri stars. 



3.3.1. Magnetic Coupling to the Disk and the Connection 

State 

The effect of the opening of the magnetic field lines 
depends on the strength of the magnetic coupling to the 
disk and how strongly the field lines can be twisted un- 
til they are severed. The variable 7(7*) = B^/B^ de- 
scribes the twisting of the magnetic field between the 
star and disk. This twisting occurs rapidly so a steady 
state configuration depends on how well the field cou- 
ples to the disk (i.e., the balance between the differential 
rotatio n and the tendency for the magnetic field to un- 
twist). lUzdenskv et al.l ([2002D describe this coupling by 
a dimensionless magnetic diffusivity parameter, 



/3 



Vt 



(17) 



where rjt is the effective magnetic diffusivity and is of 
the order of magnitu de of the disk's effective viscosity 
(|Lovelace et al.|[l995l ). H is the scale height of the disk, 
and Vk is the Keplerian rotation velocity. MP05 assume 
/3 is constant throughout the disk. The field is strongly 
coupled to the d isk for values of / 3 < 1 and weakly cou- 
pled for (3 > 1. lUzdenskv etafl (|2002t ) find that when 
7 exceeds a value of order unity (defined by the critical 
twist parameter j c ) the magnetic field will be severed be- 
cause the magnetic pressure force associated with B$ will 
push outward and cause the dipole field loops to open. 
The magnetic field is connected to the disk only in the lo- 
cation where I7I < j c . MP05 use the values j3 = 0.01 and 
7 C = 1 in their models, and we adopt the same fiducial 
values in this work. They suggest that j3 = 0.01 is the 
most probable value for a T Tauri accretion disk with the 
use of an a model prescription; however it is uncertain 
that disks surrounding massive stars will have this same 
value. For example, massive stars emit more ionizing ra- 
diation which will yield a higher ionization fraction on 
the disk surface, causing j3 to decrease, but these disks 
are also more massive than those surrounding low-mass 
PMS stars and are therefore thicker, causing j3 to in- 
crease. To account for our uncertainty in this parameter 
we experiment with different values in the following sec- 
tion. 

MP05 show that the magnetic connection between the 
star and disk changes at a threshold value of the stellar 
spin rate. Specifically, the stellar magnetic field will only 
be connected to a small region of the disk within i? co if 
the stellar rotation rate as a fraction of break up, 
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Fig. 3. — The top left panel shows the stellar radius as a function 
of stellar mass for stars with masses 0.5 - 50 Mq. The other panels 
show the stellar radius (top right), stellar period (bottom left), 
and stellar spin rate as a fraction of break up (bottom right) as a 
function of time for stars with masses 0.5 - 50 Mq. Figure [2] shows 
the accretion histories. 



falls below: 



where 



/< 1 



Pic) (lc~>P) , 



2 p5/2 



(19) 



(20) 



is a dimensionless parameter that relates the strength of 
the magnetic field to the accretion rate. This connection 
state, which MP05 denote as state 1, will result in no 
spin-down torques transferred to the star. If / exceeds 
this value then the system is in state 2 which is charac- 
terized by a magnetic connection on either side of R co 
resulting in both spin-up and spin-down torques acting 
on the star. 



-7/2 



1 - 



Rt_ 

Rco 



3/2' 



P 



(22) 



We assume the accreted disk material is quickly inte- 
grated into the structure of the star and adds angular 
momentum to the star at a rate given by equation (6) 
where Ra is replaced by Rt- This material acts to spin 
up the star. 

The magnetic connection over a range in radii in the 
disk can extract angular momentum from the star and 
transfer it to the disk. If the system is in state 2 then the 
magnetic field is connected to the disk from R t to i? out = 
(1 + /?7c) 2 ^ 3 R co which yields a net magnetic torque on 
the star: 



2 (1 + /37c)" 1 - (1 + Plc)~ 



-2(R co /R t f 2 + (R co /R t f 



(23) 



If the system is in state 1 then the magnetic field is con- 
nected to only a small portion of the disk which leads to 
a n egligible torque on the star, so we set r m = follow- 
ing [MatTeTaD ((2010f ). Note that equation (23) reduces 
to equation (5) for the limiting case of no field opening 
(7 C — > oo), marginal coupling (/3 = 1), and a disk that 
is truncated at the Alfven radius (Ra) and extends to 
infinity. 



Parameter 


Fiducial Value 


£ 


1 g cm~ 2 


M D 


0.02 M*, f 


t„ 


10 6 yr 


B* 


2 kG 





0.01 


7c 


1 



TABLE 1 

Table of fiducial values used for our model parameters. 



3.3.2. Magnetic and Accretion Torques 

The twisting of the magnetic field by the differential 
rotation between the star and disk causes torques to be 
conveyed between the two. The twisting of the mag- 
netic field within R co leads to spin-up torques whereas 
the field lines connected to the disk outside of i? co act to 
spin down the star. If the magnetic field is strong enough 
then the disk will be disrupted by the stellar magneto- 
sphere where the magnetic stress is able to maintain the 
accretion rate within the disk. At this location, denoted 
by Rt , the magnetic stress is large enough to remove the 
excess angular momentum and funnel the disk material 
along the magnetic field lines. This material and its an- 
gular momentum is transferred to the star. If R t > R co 
the magnetic stress hinders the accretion rate. 

The location of Rt depends on the connection state of 
the system. In state 1 the truncation radius is 

Rt = ( 7 c^) 2/7 i?*. (21) 
In state 2 the truncation radius is given by 



4. RESULTS 

The initial star-forming core properties are determined 
by the core mass (M core ), core density profile (k p ), and 
core surface density (S). These parameters control the 
accretion rate for the primary accretion phase as de- 
scribed in §3.2.1. We initially create a "pre-collapse" ob- 
ject with a mass less than 0.01 Mq which grows in mass 
with the accretion rate given by equation (13). When the 
object reaches a mass of 0.01 Mq we initialize our pro- 
tostcllar and angular momentum evolution model and 
assume the protostar is initially rotating at 1% of its 
break up speed. When the protostar is initialized, it is 
immediately spun up since the accretion rate is large, so 
our chosen value for the initial rotation speed is unim- 
portant. We solve equa tion (11) wi t h the fourth-order 
Runge-Kutta scheme of iPress et alJ (|2007t ) and update 
the angular momentum of the star by computing the 
net torque on the star arising from the accretion and 
magnetic torques described in §3.3.2. We use this re- 
sult to update Q+. We cap the stellar rotation rate at 
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50% of breakup, a limit imposed by gravitational torques 
(|Lin et al.ll2011l) . The fiducial values used for our model 
parameters are given in table [TJ 

4.1. Effect of the Star-Disk Magnetic Interaction 

Figure [3] shows the radial and rotational evolution for 
stars ranging in final stellar mass from 0.5 - 50 Mq. 
These models were simulated with the fiducial param- 
eters given in table [T] The disk-clearing accretion phase 
is assumed to last 3 Myr, although as discussed in § 
1, this assumption is almost certainly not correct for 
high mass stars. As we show below, using a shorter 
disk clearing timcscale for the massive stars would only 
strengthen our results. We choose to run the disk clear- 
ing phase for three decay time scales because accretion 
disks around low-mass st ars survive for several million 
years (jHerbst et al.l 120071 ) , with an accretion rate that 
likely decreases with time. The swelling in radius by a 
factor of three, shown in the upper plots of Figure [31 
is a result of t he star transitioning from a convective to 
radiative core (Hosoka wa fc Omu kai 2009), which redis- 
tributes entropy within the star. For the stars presented 
in Figure [3J this occurs in the primary accretion phase 
for the most massive stars {M+j > 15 M Q ) and during 
the disk clearing accretion phase for the 0.5, 1, and 5 M Q 
stars. If the jump in radius occurs during the main accre- 
tion phase, it causes the star to immediately slow down, 
but the star is almost instantly spun back up because of 
the high accretion rate. In the case of the 5 M Q star, 
this jump in radius also significantly decreases the spin 
rate of the star, but since it occurs when the accretion 
rate is much lower the star only gradually spins up as it 
contracts and accretes material. In contrast, for the 0.5 
and 1 Mq stars magnetic torques are able to continue to 
spin dow n the star a f ter th e jump in radius occurs. We 
note that lMatt et alJ (|2010l ) produced Sun-like stars with 
faster rotation rates (~ 20 — 40% of break up) perform- 
ing a similar analysis. We report a lower rotation rate 
for our 1 M Q protostar because it h as a differen t radia l 
history than the stars produced bv iMatt et alJ (|2010f ). 
Our 1 Mq protostar contract s more slowly tha n the 1 
M protostar model used bv IMatt et"aT1 l(2010f ). After 
3 Myr, our model gi ves a radius of 3.8 R© as compared 
to IMatt et all l|2010D 's ~ 3 R . At times <1 Myr, the 
model radii can differ by factors of ~2. The larger radii 
in our model produce more spin-down. The differences 
in predicted radii likely arise because our model accounts 
for the extra entropy provided both by deuterium bur n- 
ing and by ongoing accretion, while IMatt et al.l (|2010D 's 
does not. We do warn, however, that there are signifi- 
cant uncertainties in how much of the accretion entropy 
is actually absorbed by the star, and differing assump- 
tions on this poi nt can produce signific ant differences in 
radial evolution (jHosokawa et al.ll2011f ). 

We find that the torques that arise from the star-disk 
magnetic interaction are unable to spin down both low- 
mass and massive protostars during the main accretion 
phase, but are important during the disk clearing phase, 
especially for low-mass stars. Low-mass stars begin to 
spin down the instant the disk clearing accretion phase 
begins whereas it takes approximately 2 Myr to begin 
to spin down massive stars for our chosen fiducial values. 
This suggests that massive stars are difficult to spin down 
due to their larger inertia and because their magnetic 
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Fig. 4. — Snapshots of the stellar radius (upper left), disk trun- 
cation radius (upper right), stellar period (lower left), and rotation 
rate as a fraction of break up (lower right) as a function of stellar 
mass taken at different times during the disk clearing phase for our 
fiducial case. The times in the legend represent the time that has 
elapsed since the disk clearing phase began. The points located in 
the bottom right panel represents the minimum mass of stars ro- 
tating at > 20% of its break up speed. We use this as an indicator 
of the transition between slow and fast rotators. 



fields are weaker relative to their stellar binding energy 
as compared to low mass stars. 

Figure 2] shows snapshots of the stellar radius, disk 
truncation radius, stellar period, and stellar rotation rate 
as a fraction of break up as a function of stellar mass 
taken at different times during the disk clearing phase. 
First consider the upper left panel, showing radius ver- 
sus mass at different times. The R— M relation toward 
which the models converge at high mass is the ZAMS; by 
3 Myr all stars above ~ 2 Mq have reached it. At smaller 
masses, the maximum radius occurs at a mass that cor- 
responds to stars that have just made the convective- 
radiative core transition at a given time. This value shifts 
to progressively smaller masses at later times. 

An interesting feature of Figure [4] is that the stellar ro- 
tation rates as a fraction of break up show a bimodal dis- 
tribution: stars with M* f < 1M© rotate at ~10% of their 
break up speed whereas stars with M*.f > 6M Q are rapid 
rotators. In between these plateaus (i.e., the "transition 
region" ) the rotation rates as a fraction of break up in- 
creases with stellar mass. Furthermore, as time increases 
we find that the ratio of rotation speed to break up speed 
decreases on both plateaus, but that this decrease is more 
noticeable for the fast rotator plateau. This is because 
the stars located on the fast rotator plateau have already 
reached the ZAMS and are no longer contracting whereas 
those located on the slow rotator plateau are easy to spin 
down because of their low inertia, even though they are 
still contracting towards the ZAMS. In contrast, we find 
that the rotation rates as a fraction of break up of the 
stars in the transition region increases with time. This 
suggests that the magnetic torques conveyed by the star- 
disk interaction are unable to counteract the increase in 
the stellar spin rate due to contraction for stars in the 
transition region. However, once these stars have reached 
the ZAMS magnetic torques do become important. The 
points located in the bottom right panel of Figure [4] rep- 
resents the minimum mass of stars rotating at >20% of 
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Fig. 5. — Same as figure 4 but all quantities arc shown at a time 
of 3 Myr, and we vary S as indicated in the legend. 

their break up speed. We use this as an indicator of 
the transition between slow and fast rotators, which we 
discuss further in §4.3. 

4.2. Sensitivity to Model Parameters 

In the previous subsection we found that massive stars 
are much more difficult to spin down than low-mass stars. 
This causes low-mass stars to become slow rotators and 
massive stars to be rapid rotators, yielding a bimodal 
distribution in stellar rotation speeds as a fraction of the 
break up speed. To explore if this qualitative result is 
sensitive to our chosen model parameters, we vary certain 
parameters while holding the other parameters fixed. In 
the figures that follow we see that by varying certain 
parameters we do not lose this feature, but only alter it. 



4.2.1. 



ng E 



Figure [5] shows the final stellar radius, disk truncation 
radius, stellar period, and rotation rate as a fraction of 
break up as a function of final stellar mass for different 
values of the initial core surface density, E. The accretion 
rate during the main accretion phase increases for higher 
£, so varying this value affects the accretion history only 
during this phase. We find that this parameter has little 
to no effect on the final spin rate of the stars because the 
magnetic torques are unimportant during this accretion 
phase. The very minor differences that do appear arise 
because the value of £ affects the time at which a star of a 
given final mass reaches the swelling phase: the swelling 
phase of the star occurs earlier in time at lower E. For 
each value of E used in our models there is a slight kink 
in between w 3 — 6 M Q and the location of this 

kink decreases in mass for smaller values of E. Stars 
to the right of this kink experience the jump in radius, 
discussed in §4.1, before the end of the main accretion 
phase whereas those to the left experience the swelling 
during the disk-clearing accretion phase. However, the 
net effect on the stellar rotation rate is obviously minor. 

4.2.2. Varying M D 

Figure [5] shows the final stellar radius, disk trunca- 
tion radius, stellar period, and rotation rate as a fraction 
of break up as a function of final stellar mass for dif- 
ferent values of the initial disk mass, Md, used for the 
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Fig. 6. — Same as figure 4 but all quantities arc shown at a time 
of 3 Myr, and we vary Md as indicated in the legend. 
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Fig. 7. — Same as figure 4 but all quantities arc shown at a time 
of 3 Myr, and we vary t a as indicated in the legend. 

disk clearing accretion phase. Increasing Md increases 
the accretion rate during the disk clearing phase, thus 
increasing the accretion torque. A larger accretion rate 
also causes the disk to be truncated closer to the star, 
effectively reducing the net spin down magnetic torque. 
This is because the stellar magnetic field lines will con- 
nect to a greater portion of the disk within R co yielding 
greater spin up magnetic torques on the star while the 
magnetic spin down torques remain unchanged. We find 
that altering Md changes the location and shape of the 
transition between the slow and fast rotation plateaus, 
but the qualitative result that rotation rates are bimodal, 
with slow rotation at low mass and rapid rotation at high 
mass, remains unchanged. Also note that the models 
converge in the limit Md — > 0. 

4.2.3. Varying t a 

Figure [7] shows the final stellar radius, disk truncation 
radius, stellar period, and rotation rate as a fraction of 
break up as a function of final stellar mass for different 
values of the disk decay time scale, t a , used for equation 
(16) . Smaller values of i a , as compared to our fiducial 
value of 1 Myr, correspond to a higher initial accretion 
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Fig. 8. — Same as figure 4 but all quantities arc shown at a time 
of 3 Myr, and we vary as indicated in the legend. 

rate that declines more rapidly for the disk clearing ac- 
cretion phase. This yields lower final spin rates at the 
end of 3 Myr. However, the overall shape of the distribu- 
tion of final spin rates as a function of stellar mass does 
not change. 

4.2.4. Varying 

Figure [8] shows the final stellar radius, disk truncation 
radius, stellar period, and rotation rate as a fraction of 
break up as a function of final stellar mass for differ- 
ent values of the stellar magnetic field strength. Clearly, 
a larger magnetic field strength provides a greater spin 
down torque on the star, yielding smaller final spin rates 
as a function of mass. As can be seen in this figure, stars 
above M+j > 2 M Q require surface fields greater than 
1 kG to experience any significant spin down torques 
and do not become slow rotators, f^/^bu ^ 0.1, un- 
less the field reaches ^10 kG. Magnetic fields this large 
have only been detected in th e chemically peculiar (e.g. , 
helium strong) Ap /Bp stars (jBorra fc Landstreetl 119791 : 
lOksala et al.ll2010() . Generally, as the field strength in- 
creases the final spin rates decrease, but the qualitative 
division between slow and fast rotators remains. We also 
find that this same trend in rotation rates as a fraction of 
break up occurs as the field lines become weakly coupled 
to the accretion disk, while holding the magnetic field 
strength fixed, as discussed next. 

4.2.5. Varying j3 and 7 C 

Figure IHl shows the final stellar radius, disk truncation 
radius, stellar period, and rotation rate as a fraction of 
break up as a function of final stellar mass for different 
values of j3 and 7 C . These parameters describe the cou- 
pling and connection of the stellar magnetic field lines 
to the accretion disk (i.e., the location where the field 
lines open and disconnect from the disk). A larger j3 
for a given 7 C increases the extent of the connected disk 
region. This is because the coupling of the stellar field 
lines to the disk acts to resist the twisting of these lines 
due to the differential rotation between the star and disk. 
Thus, weaker field coupling will lead to a greater spin- 
down torque acting on the star leading to lower rotation 
rates as depicted in Figure [9j Likewise, a greater j c for 
a given f3 will allow the field lines to experience a greater 



Fig. 9. — Same as figure 4 but all quantities arc shown at a time 
of 3 Myr, and we vary /3 and f c as indicated in the legend. 

twist before opening, also increasing the size of the con- 
nected disk region. For the case where 7 C — » oo (i.e., 
field lines are allowed to twist to large values without 
opening), the field lines will connect to the whole disk 
outside i?t- This will lead to a greater spin down torque. 
The case where j3 = 1 and j c = oo reduces to the case 
described in section 2. Figure [9] shows that as (3 increases 
for 7 C = 1, all stars have lower rotation rates. However 
the two plateaus still remain. 

4.3. The Characteristic Mass for the Slow to Fast 
Rotator Transition 

In this work we have found a robust division between 
slow and fast rotators. Specifically, we find that low-mass 
stars (e.g., stars with M* < 1 M Q ) are slow rotators, 
easily spun down via magnetic torques that arise from 
the star-disk interaction, and rotate at ~ 10% of their 
break up speed, whereas massive stars (e.g., M+ > 6 M ) 
are preferentially fast rotators. This is because massive 
stars are difficult to spin down due to their larger iner- 
tia and because their magnetic fields are weaker relative 
to their stellar binding energy as compared to low mass 
stars. Furthermore, this division is also dependent on 
the R-M relationship. The stars located on the fast ro- 
tator plateau have reached the ZAMS by the end of the 
main accretion phase or early on during the disk clear- 
ing phase; whereas, the stars located on the slow rotator 
plateau are shrinking towards the ZAMS for the entirety 
of the disk clearing phase. Likewise, the stars located in 
the transition region are contracting towards the ZAMS 
for a significant portion of the disk clearing phase but are 
contracting much faster than the low-mass slow rotators, 
leading to the sudden rise in rotation rates as a fraction 
of break up. 

To further illustrate the division between slow and fast 
rotators for each of our model parameters, in figure [TU] 
we plot the minimum stellar mass at which the star ends 
accretion rotating at 20% of its break up speed, which we 
call A/20- Each panel shows how A/20 depends on the in- 
dividual parameters in our model (while setting the other 
parameters to their fiducial values). The top panels show 
that M 2 o decreases by only a small amount as the disk 
lifetime (i.e., the amount of time the disk survives and 
supplies mass to the star during the disk clearing phase) 
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or disk decay time scale increases. We also see that this 
characteristic mass, as a function of the initial core sur- 
face density, is relatively constant as indicated by the 
nearly horizontal line on the middle right panel of fig- 
urc llOl In contrast, M20 spans a larger mass range as we 
vary the initial disk mass used for the secondary accretion 
phase as shown in the middle left panel. This is because 
the accretion rate, and therefore the accretion torque, is 
proportional to the disk mass used in our model. We 
find that as Md — > the values of M20 become constant 
but we notice that M20 decreases most as the disk mass 
increases from ~ 10~ 3 — 10~ 2 M*. The division between 
slow and fast rotators slowly decreases in stellar mass 
for initial disk masses above ~ 10 -2 M*. Even though 
varying this parameter leads to larger variations in M20 
as compared to the top panels, it does not change the 
qualitative division between slow and fast rotators. 

The bottom panels of figure [10] show how the slow-fast 
rotator division is affected by the stellar magnetic field 
strength and the coupling of the stellar magnetic field 
lines to the disk, which are the parameters that are re- 
sponsible for the removal of angular momentum from the 
star. The black solid lines in these panels show that the 
division between the slow and fast rotators (i.e., M20) 
diverges for large magnetic field strengths > 4.5 kG) 
or weak held coupling (/3 > 0.05) for a disk clearing ac- 
cretion phase that lasts for 3 Myr. This is because no 
stars will be rotating at or above 20% of their break up 
speed at the end of 3 Myr for such high values of 
or p. For comparison, and also because we expect disks 
to have shorter lifetimes around massive stars, we also 
include the values of M 2 o at 0.5 Myr after the disk clear- 
ing phase began (teal dotted lines). We find that M 2 o 
is larger at shorter times because these stars are still 
contracting towards the ZAMS. At 0.5 Myr stars with 
masses greater than ~ 5 M Q have reached the ZAMS, as 
indicated by the kink and faster increase of M20 in these 
plots for the 0.5 Myr case. 

5. DISCUSSION 

We have shown that massive stars are fast rotators at 
birth and that their initial rotation rates are unlikely 
to be regulated by the star-disk magnetic interaction. 
We have found that magnetic torques can only effec- 
tively spin down massive stars that have low accretion 
rates, long disk lifetimes, weak magnetic coupling with 
the disk, and/or surface magnetic fields that are signifi- 
cantly larger than what current observational estimates 
suggest. We thus conclude that their initial rotation 
rates are likely regulated by gravitational torques. Since 
massive stars arrive on the main sequence as fast rota- 
tors, their variation in rotation rates as a fraction of their 
break-up rate is likely a result of evolutionary spin down, 
due to stellar expansion and/or angular momentum loss 
via stellar winds while on the main sequence. 

5.1. Observational Implications 

A topic of current debate is whether the distribu- 
tion of the projected rotational velocities of massive 
stars depend on birth environment or if t his property 
is on ly affected by evolutionary spin down ( Strom et al.l 
20051: iDufton et al.l [20061: iHuang fc Giesl 20061 120081: 
Wolff et al.ll2007l 120081: IHuang et al.ll2010D . Istrom et all 
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Fig. 10. — This figure illustrates the sensitivity of the model 
parameters. The y-axes show the minimum stellar mass where 
/ > 0.2, denoted as M20, for different parameters as indicated on 
the x-axes. Except for the top left plot, the black solid lines indicate 
that these values were taken for a disk lifetime of 3 Myrs. In the 
bottom panels, the teal dashed lines show the value of M20 0.5 Myrs 
after the beginning of the disk clearing phase for comparison. 



(j2005| ) observed the rotational velocities of B stars lo- 
cated in high stellar density clusters and compared them 
to field stars of similar age (~12-15 Myr). They found 
that, on average, the cluster stars had larger rotational 
velocities than the field stars in their sample and that 
only the most evolved cluster stars had similar rota- 
tional VBtocitiesas their field star counterparts. Like- 
wise, iWolff et all ([20071 12008ft observed that massive 
stars (e.g., M* > 6 Mq) found in clusters characterized 
by a high stellar density are faster rotators than their 
similar mass counterparts located in lower density clus- 
ters. These studies concluded that the initial spin rates 
of these stars depend on the initial star-forming envi- 
ronment since these stellar ensembles, which have sur- 
vived as bound clusters, likely form in molecular clouds 
characterized by high surface densities. Furthermore, 
IWolff et al.l (|2007t ) compared the distribution of the rota- 
tional velocities of B stars in both young and older high 
density and low density environments and did not detect 
a significant evo lutionary change. 

In agreement. IHuang et al.l ([2010f ) compared the rota- 
tion rates of cluster and field B stars and found that, 
on average, cluster stars tend to rotate faster than field 
stars. However, by grouping the stars by surface gravity, 
an age proxy, they found there is little difference between 
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the average rotational velocities for the field and cluster 
stars as a function of age, and that they exhibit a simi- 
lar spin-down with advanced evolution. They also found 
that field stars are in general more evolved than cluster 
stars. These results suggest that the observed trend in 
the rotational velocities of B stars are due to evolution- 
ary spin down rather than to the initial conditions of 
the environment in which they formed. They argue that 
the discrepancy between the average rotation rate of the 
field stars and cluster stars in their sample is that the 
field stars have undergone evolutionary spin down since 
the field star sample contained more evolved stars. 

For a fixed surface magnetic field strength, we find here 
that the initial rotation rates of massive stars, due to disk 
locking, have no dependence on the environmental den- 
sity. As described in §3.2.1, the accretion rate during 
the main accretion phase does depend on the star form- 
ing environment, with larger surface density yielding a 
greater time-averaged accretion rate. iWolff et al.l (|2007f ) 
proposed that the higher rotation rates they report for 
stars in dense clusters are the result of disk-locking plus 
a systematically higher accretion rate in dense clusters. 
However, we find that magnetic torques are insignificant 
during the main accretion phase regardless of environ- 
ment density because of the high accretion rates. These 
torques only become important during the disk clearing 
phase, and there is no obvious reason that the properties 
or behavior of the disk during this phase should depend 
on the environment. However, this does not rule out 
other factors that may depend on the environment. In 
this work we assumed that all stars had the same surface 
magnetic field strength. If the strength of the magnetic 
fields present during the star formation process depends 
on environment, either because the star- forming cloud 
has a different magnetic mass to flux ratio and/or be- 
cause the ambipolar diffusion process depends on den- 
sity, then this could provide a viable explanation for the 
difference in rotational velocities of young stars in en- 
vironments of varying density. Another possibility for 
the difference in rotational velocities of stars born in dif- 
ferent environments may be related to the lifetimes of 
disks in such environments. We have found that the ro- 
tation rates of these stars depend crucially on the life- 
time of the accretion disk. Thus, if disks have shorter 
lifetimes in higher stellar density environments, possibly 
due to tidal dissipation from interactions with neighbors 
or rapid photoevaporation due to radiation from nearby 
massive stars, then the initial rotation rates of these stars 
will only increase as they contract towards the ZAMS 
(|Wolff et al.ll2007D . 

5.2. Future Work and Caveats 

In this work, we have omitted two potentially impor- 
tant effects: that magnetic fields might be stronger early 
in stars' lives, and that stars can be spun down by winds 
on the main sequence. As mentioned in §1, magnetic 
fields in massive stars are likely to be the decaying rem- 



nants of magnetic flux swept up during the star formation 
process. Therefore, it is plausible that accreting massive 
stars have stronger magnetic fields than those we ob- 
serve as main sequence O and B stars. If this is the case, 
then massive stars will likely be spun-down via magnetic 
torques. If the decay process is the same for all stars then 
we expect that the strongest magnetic fields should be 
observed in the slowest rotators. However, we also dis- 
covered that the spin rates of these stars depend heavily 
on how well the stellar magnetic field lines couple to the 
accretion disk. As described in §3.3.1 the true value of 
/3 is highly uncertain because it depends on the micro- 
physics of the accretion disk. Since observations of disks 
around massive stars are rare we are unable to provide a 
confident estimate for f3. However, by exploring a range 
of values for f3 we have determined that if the field lines 
arc weakly coupled to the disk then magnetic torques 
can sufficiently spin down massive stars. Also, measur- 
ing the rotation rates of young, massive stars can provide 
a better estimate for /?. If the slowest rotators prove to 
have weak magnetic fields, then it may be likely that the 
field lines were weakly coupled to the disk, resulting in a 
larger /3, thus producing these slower rotators. 

Stars on the main sequence also shed mass and an- 
gular momentum via stellar winds, which we have ne- 
glected in this work. In the presence of a stellar mag- 
netic field, these winds will couple with the field lines 
causing the star to lose a significant amount of angular 
momentum as it evolves ijWeber fc Davis! Il967f ). Since 
the mass loss rates of stars increa ses with stellar mass 
(|Nieuwenhuiizen fc de Jagerl Il990f ) more massive stars 
will lose angular momentum at a greater rate. If the 
spin rates of massive stars are regulated by gravitational 
torques rather than magnetic torques produced by the 
star-disk magnetic interaction then we expect that all 
massive stars should be rotating at ~50% of their break 
up speed once they are deposited on the ZAMS, assum- 
ing that their disks survive long enough. Spin down 
will occur as they evolve and shed angular momentum 

vi a stellar wind s. This is consistent with the results 

of iHuang et alJ (|2010f ) who found that young stars with 
masses greater than ~ 2 M Q are preferentially fast ro- 
tators and that the average rotation speed as a fraction 
of the break up speed, for each mass bin, decreases for 
increasing stellar mass. 
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